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Abstract 

We study the physical properties of MnAs under strain by using 
accurate first-principles pseudopotential calculations. Our results pro- 
vide new insight on the physics of strained multilayer that are grown 
epitaxially on different lattice mismatched substrates and which are 
presently of interest for spintronic applications. We compute the strain 
dependence of the structural parameters, electronic bands, density of 
states and magnetization. In the region of strain/stress that is eas- 
ily directly accessible to measurements, the effects on these physical 
quantities are linear. We also address the case of uniaxial stress induc- 
ing sizeable and strongly non linear effects on electronic and magnetic 
properties. 

Keywors: MnAs, strain, pseudopotential, magnetic properties, elec- 
tronic properties, epitaxy 

PACS: 71., 71.15.Hx, 71.15.Mb, 71.20.-b, 75.50.-y 
1. Introduction 

The recent discovery of giant magneto-resistance has driven increasing in- 
terest to the study of magnetic materials. The main target is the design 

°* Corresponding Author: Tel. +39 040 2240242; Fax. +39 040 224601; E-mail: 
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and the realization of new devices with materials whose magnetic properties 
could be tailored in accordance to the needs of electronic industry. Within the 
solid state community is rapidly becoming popular the word "spintronic" to 
denote electronic-like heterostructures where the relevant physical quantity 
is the spin of the carriers and its interactions with external magnetic fields 
rather than the charge of holes and electrons and the associated electronic 
properties. In this context the most promising materials are Mn compounds. 

Although they naturally occur in the hexagonal structure whose symme- 
try is completely different with respect to the zincblende structure of conven- 
tional semiconductors, high-quality epitaxial Mn compounds can be grown 
on many semiconductor substrates, different for composition and orienta- 
tion. Realistic cases include the growth on (lll)-As-terminated GaAs or 
(111) Si [p], where the MnAs c-axis is parallel to the (111) zincblende sub- 
strate axis, and on (001) GaAs or Si, with two possible different orientations 
(usually referred as type-A and type-B) of the MnAs c-axis with respect 
to the zig-zag anion-cation bond chain underlying the substrate surface. A 
recent experiment || has proved the possibility of growing ordered films 
(heterostructures) of MnAs on ZnSe. The presence of a ZnSe buffer layer 
exclusively stabilizes a particular phase (type-B) of MnAs; at variance, the 
coexistence of different phases in the magnetic epitaxial layers is experimen- 
tally observed on GaAs substrates 

Due to the lattice mismatch with the substrate, which is different accord- 
ing to its composition and orientation, we expect the Mn-compound epilayers 
to grow in different kind of strained structures, some of them with consider- 
able deformation with respect to the original free-standing structure, and in 
general of reduced symmetry. 

After a description of the structural, electronic and magnetic properties 
of MnAs in the free-standing hexagonal and zincblende structures, in this 
paper we address in particular the effects on these properties induced by a 
uniaxial stress along the c-axis of the hexagonal structure. 

2. Computational details 

Our calculations are performed within the spin-polarized density functional 
theory using ab-initio pseudopotentials and plane-wave formulation. For the 
exchange and correlation energy we use the generalized gradient approxi- 
mation (GGA) with the Perdew-Wang parameterization [f|], that has been 
shown to provide a good description of the metallic bond. From the com- 
parative tests that we have performed using the local density approximation 
(LDA) with the Perdew-Zunger parametrization of the exchange and correla- 
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tion functional || , we found that GGA slightly improves the agreement of our 
results with respect to the available experimental data (see Tab. 1). We have 
generated separable pseudopotentials in the Kleinman-Bylander form ||: the 
As pseudopotential has been generated according to the Martins- Troullier 
technique J7|, while the Mn pseudopotential according to the Vanderbilt 
method || which allows to have reliable results with a reasonable size of 
the plane-wave basis set (corresponding in the present case to a kinetic en- 
ergy cutoff of 30 Ry). Reciprocal space integration in the Brillouin zone is 
performed using special k-point techniques [§, 10|. A particular care in the 



choice of the sampling k-point set must be reserved to the system studied: at 
zero stress (for the free-standing hexagonal structure) we use a (12 x 12 x 8) 
Monkhorst-Pack [TO] grid that corresponds to 336 k-points in the irreducible 
wedge of the first Brillouin zone, and ensures an accurate convergence of 
all the electronic properties studied. For the strained structures this set is 
modified accordingly. The integration over the points of the Brillouin zone 
are performed with smearing techniques [ll|] using a Gaussian broadening of 

0. 02 Ry. 

Deformations can be either introduced by constraining the structures in 
some directions and allowing them to relax in the others (i.e., fixing some 
components of the strain tensor e^, as it is the case for epitaxial structures 
that must accommodate the lattice mismatch with the substrate), or by 
applying a stress cr^ (external force). We focus here on the last case, sim- 
ulating the application of a uniaxial stress along the c axis of the hexagonal 
structure, here referred to as the z direction. The correspondent equilibrium 
lattice structures must have vanishing stress along the other directions, i.e. 
°~zz 7^ and a xx = o yy = 0. In practice in the numerical experiments we con- 
sider several structures with different strain tensors e^, with e xx = e yy ^ e 22 , 
and we compute the stress tensor according to the formalism introduced 
by Nielsen and Martin. [12, |13| The structures can be considered optimized 
if the residual stress in the xy plane is within RilkBar. We point out that 
in general this condition is not equivalent to the one of volume conservation, 

1. e. zero-trace strain tensor, and requires a larger computational effort. 



3. Results 

3.1 Free-standing unstrained structures 

In normal free-standing conditions MnAs is a ferromagnetic metal with the 
NiAs structure, i.e. hexagonal lattice with four atoms in the unit cell, which 
is described by the parameters a (corner of the hexagon) and c (Fig. 1). The 
zincblende structure is experimentally observed only for the related alloy 
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Gai-zMn^As with the extreme diluted concentration of Mn atoms, x < 0.07. 
We will address here also the study of this structure, for a comparison with 
the hexagonal one. In Tab. 1 we summarize the results for the structural 
and magnetic properties obtained with our pseudopotentials with GGA, to- 
gether with those obtained within LDA, and also by other authors using 
full-potential linearized augmented plane- wave techniques (FLAPW) fl4|| . 
The comparison with the available experimental data completes the Table. 

Comparing zincblende and hexagonal MnAs, we found that the last one 
is more dense and more stable (by ~0.7 eV); the increase in the atomic 
distance in the zincblende structure is accompanied by an increase in the 
total magnetization M to t, which goes from 2.80 \xb to 3.84 /ie, to be com- 
pared with the limiting value of 4 \xb per MnAs pair for atoms at infi- 
nite distance. We compute both the total and the absolute magnetization, 



dr. 



M ^ = ft, / (?( r ) - p 1 ^) dr and M <*» = nti / k T oo - P 

The two structures are very different as far as the electronic properties 
is concerned: at variance with the hexagonal phase, the zincblende MnAs 
is half-metallic, as we show in Fig. 2. From the density of states (DOS) 
for the majority and minority spin projected onto the different atomic states 
(not shown here) it is evident the atomic origin (d states of Mn) of the flat 
bands in the two spin orientations. Our results are in agreement with some 
obtained from previous pseudopotential calculations. |I5| 



3.2 Strained structures for applied uniaxial stress 

The uniaxial stress along the c axis does not cause any reduction of the 
symmetry of the structure, which remains hexagonal but with a c/a ratio 
different from the free-standing one. For uniaxial stresses up to about 5.2 
GPa the effects on the structural parameters, electronic structure and mag- 
netic properties are almost linear. For the electronic structure we report the 
bands for the unstrained and for one strained configuration (Fig. 3, solid and 
dashed lines), and the total DOS projected onto the different atomic states 
for the same configurations and for a third one with intermediate strain con- 
figuration (Fig. 4). No dramatic changes are observed, only a progressively 
filling of the DOS around the Fermi energy. The total magnetization reduces 
linearly with the stress from 2.80 /xg to 2.65 /Jb- 

Sizeable changes occur instead for higher uniaxial stress. At about 7.6 
GPa the effect of the applied stress on the structural properties is no longer 
linear: it can be seen in terms of induced strain, or volume change (which 
amounts to —10%), or reduction of c/a ratio (of about —20%, thus corre- 
sponding to a strong cell deformation). Also the effects on the electronic and 
magnetic properties are important. The band structure (Fig. 5) looks quite 
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different from the free-standing and from the low-stress region: the bands 
originating from the d-states of Mn approach the Fermi level, going towards 
higher energy for the majority spin, and towards lower energy for the minor- 
ity spin. The net effect is a dramatic reduction of the magnetization, which 
goes to 0.26 hb- Although further investigation is needed to better under- 
stand the microscopic origin of such effect, our preliminary results seem to 
indicate that the magnetization reduction is related not simply to a sizeable 
compression, but, more important, to a strong deformation of the structure. 



4. Summary 

We calculate the variation of the structural parameters, electronic bands, 
density of states and magnetization of MnAs under the application of a uni- 
axial stress along the c axis of the hexagonal structure. We found a region 
of strain/stress where the effects are linear and rather small, but also a re- 
gion of sizeable and strongly non linear effects on electronic and magnetic 
properties. 
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Table 1: Structural parameters and magnetization for MnAs in the FCC 
and HEX phases compared with other calculations and with available exper- 
iments. The pseudopotential results, both with LDA and GGA, are those 
obtained in the present work. 
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Figure 1: The NiAs structure. 



Figure 2: Majority and minority bands for the half-metallic zincblende 
MnAs. The energy zero is set to the Fermi energy. 



Figure 3: Majority and minority bands for the ferromagnetic MnAs in 
the unstrained structure and in the structure corresponding to an applied 
uniaxial stress on c of 5.2 GPa. The energy zero is set to the Fermi energy. 



Figure 4: Total density of states projected onto the different atomic states 
for the unstrained hexagonal MnAs, and for two different strain conditions 
corresponding to an applied uniaxial stress on c respectively of 2.6 GPa and 
5.2 GPa. 



Figure 5: Majority and minority bands for the ferromagnetic MnAs in 
the unstrained structure and in the structure corresponding to an applied 
uniaxial stress on c of 7.6 GPa. The zero is arbitrary set to the Fermi energy. 
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